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Abstract — Physical layer network-coding for the n-way wireless 
relaying scenario is dealt with, where each of the n user nodes Xi , 
Xi,...,X n wishes to communicate its messages to all the other 
(n— 1) nodes with the help of the relay node R. The given scheme, 
based on the denoise-and-forward scheme proposed for two-way 
relaying by Popovski et al. in QJ, employs two phases: Multiple 
Access (MA) phase and Broadcast (BC) phase with each phase 
utilizing one channel use and hence totally two channel uses. 
Physical layer network-coding using the denoise-and-forward 
scheme was done for the two-way relaying scenario in |2|, for 
three-way relaying scenario in |3|, and for four-way relaying 
scenario in 1 11]. This paper employs denoise-and-forward scheme 
for physical layer network coding of the n-way relaying scenario 
illustrating with the help of the case n = 5 not dealt with so far. 
It is observed that adaptively changing the network coding map 
used at the relay according to the channel conditions reduces the 
impact of multiple access interference which occurs at the relay 
during the MA phase. These network coding maps are chosen 
so that they satisfy a requirement called exclusive law. We show 
that when the n users transmit points from the same A/-PSK 
(M — 2 A ) constellation, every such network coding map that 
satisfies the exclusive law can be represented by a n-fold Latin 
Hyper-Cube of side M. The singular fade subspaces resulting 
from the scheme are described and enumerated for general values 
of n and M and are classified based on their removability in the 
given scenario. A network code map to be used by the relay for 
the BC phase aiming at reducing the effect of interference at the 
MA stage is obtained. 

I. Background And Preliminaries 

The two-stage protocol for physical layer network coding 
for the two-way relay channel first introduced in |4|, exploits 
the multiple access interference occurring at the relay so that 
the communication between the end nodes can be done using 
a two stage protocol. The works in (5J, (6| deal with the 
information theoretic studies for bidirectional relaying. In (2j, 
modulation schemes to be used at the nodes for uncoded 
transmission for the two-way relaying were studied. 

The work done for the relay channels with three or more 
user nodes is given in (3), ||7)-|fTT). In (7), authors have 
proposed a two stage operation for three-way relaying called 
joint network and superposition coding, in which the three 
users transmit to the relay node one-by-one in the first phase, 
and the relay node makes two superimposed XOR-ed packets 
and transmits back to the users in the BC phase. The packet 
from the node with the worst channel gain is XOR-ed with the 
other two packets. The protocol employs four channel uses, 
three for the MA phase and one for the BC phase. It is claimed 
by the authors that this scheme can be extended to more 




Fig. 1. An n-way relay channel 

than three users as well. In the work by Pischella and Ruyet 
in (8) a lattice-based coding scheme combined with power 
control, composed of alternate MA and BC phases, consisting 
of four channel uses for three-way relaying is proposed. The 
relay receives an integer linear combination of the symbols 
transmitted by the user nodes. It is stated that the scheme 
can be extended to more number of users. These two works 
essentially deal with the information theoretic aspects of multi- 
way relaying. An 'opportunistic scheduling technique' for 
physical network coding is proposed by authors Jeon et al. 
in PU| , where using a channel norm criterion and a minimum 
distance criterion, users in the MA as well as the BC phase 
are selected on the basis of instantaneous SNR. This approach 
utilizes six channel uses in case of three-way relaying and 
it is mentioned that the approach can be extended to more 
number of users. In ||9), a 'Latin square-like condition' for 
the three-way relay channel network code is proposed and 
cell swapping techniques on these Latin Cubes are suggested 
in order to improve upon these network codes. The protocol 
employs five channel uses, and the channel gains associated 
with the channels are not considered in the construction of 
this network coding map. 

We consider the n-way wireless relaying scenario shown in 
Fig. 1, where n-way data transfer takes place among the nodes 
X%, X2,---, X n with the help of the relay R assuming that the 
n nodes operate in half-duplex mode. The relaying protocol 
consists of two phases, multiple access (MA) phase, consisting 
of one channel use during which X\, X2,..., X n transmit to R; 
and broadcast (BC) phase, in which R transmits to X\, Xi,..., 
X n in a single channel use. Network Coding is employed at 
R in such a way that Xi can decode X,'s message for i,j = 
1,2, ...,n and j 7^ i, given that Xi knows its own message. 
Latin Cubes have been explored as a tool to find the network 
coding map used by the relay, depending on the channel gain 



in |3j. The throughput performance of the two stage protocol 
for three-way relaying given in [3] is better than the throughput 
performance of the 'opportunistic scheduling technique' given 
in 1 10) at high SNR, as can be observed from the plots given 
in | 3] . The work in (TTJ further extends the approach used 
in I 3| to four-way relaying and employs two channel uses 
for the entire information exchange amongst the four users, 
which makes the throughput performance of the scheme better 
than the other existing schemes. This scheme that utilizes two 
channel uses for the entire information exchange between three 
and four users using a relay in |3) and (JTT) respectively, is 
extended to n users in this paper, for the case when Af-PSK 
is used at the end nodes. 

For our physical layer network coding strategy we use the 
mathematical structure called a Latin Hyper-Cube defined as 
follows: 

Definition 1: An n-fold Latin Hyper-Cube L of r-th order 
of side M fL?) is an M x M x ... x M (n times) array 
containing M n entries, M n ~ r of each of M r kinds, such that 
each symbol occurs at most once for each value taken by each 
dimension of the hyper-cube. [|] 

For our purposes, we use only n-fold Latin Hyper-Cubes of 
side M on the symbols from the set Z t = {0, 1, 2, t — 1}, 
t > M"" 1 . 

II. Signal Model 

Multiple Access (MA) Phase: 

Suppose user node X k for k = 1,2, ...,n wants to send 
a A-bit binary tuple to all user nodes Xi for / = 1,2, ...,n 
and I ^ k. The symmetric Af-PSK (M = 2 A ) constellation, 
given by S = {e 2 * k / M \k = 0, 1, .... A/ - 1} is used at X x , 
X 2 ,..., X n , and /1 : Fj ^ 5 denotes the map from bits to 
complex symbols used at X\, X 2 ,..., X n where F2 = {0, 1}. 
Let x\ = fi (si) , x% = p (s 2 ) ,...,!„ = /i (s„) G S denote the 
complex symbols transmitted by X\, X 2 ,..., X n respectively, 
where si, s 2 , s n € ¥ 2 . Here, we assume that the Channel 
State Information (CSI) is not available at the transmitting 
nodes and perfect CSI is available at the receiving nodes. The 
received signal at R in the MA phase is given by, 

Y R = Hixi + H 2 x 2 + ... + H n x n + Z Rl (1) 
where H\, H 2 ,..., H n are the fading coefficients associated 
with the Xl-R, X 2 -R,..., X n -R link respectively. The additive 
noise Z R — CAT (0, a 2 ), where CAf(0,a 2 ) denotes the 
circularly symmetric complex Gaussian random variable with 
variance a 2 . 

The effective constellation seen at the relay during the MA 
phase, denoted by S R (Hi, H 2 , H n ), is given by, 
Sr (Hi, H2, H n ) = {H\x\ + H2X2 + ... + H n x„\xi,X2, ...,x n e S} 

The minimum distance between the points in the constel- 
lation S R (Hi,H 2 ,...,H n ) denoted by d mm (H X ,H 2 , ...,H n ) 
is given in ([3]l on the next page. From d3), it is clear 
that there exists values of (Hi, H 2 , H n ), for which 
(H\,H 2 , ...,H n ) — 0. 

'The definition has been modified slightly from the referred article "On 
Latin and Hyper-Graeco-Latin Cubes and Hyper Cubes" by K. Kishen 
(Current Science, Vol. 11, pp. 98-99, 1942), in accordance with the context. 



Definition 2: A fade state (Hi, H 2 , H n ) is defined 
to be a singular fade state for the MA phase of n- 
way relaying, if d mm (Hi, H 2 , H n ) = 0. Let U = 
{ (Hx , H 2 , H n ) G C" I d min (Hi , H 2 , H n ) = 0} denote 
the set of all singular fade states. For singular fade states, 
\S R (Hi,H 2 ,...,H n )\<M n . 

Let the Maximum Likelihood (ML) estimate of 
(xi,x 2 , ...,x n ) be denoted by (ii,x 2 , ...,x n ) € S n at 
R based on the received complex number Y R , i.e., 

(xi,&2, — ,x n ) = axg min \\Y R -HX\\, (2) 

where H — [Hi H 2 ... H n ) and X — [xi x 2 ... x n ] T . 
Broadcast (BC) Phase: 

During the BC phase, the received signals at 
Xi , X 2 , X n are respectively given by, 

Y Xk =H' k X R + Z k , fc = l,2,...,n; (8) 

where X R = M Hl - H '^-' Hn ((xi,x 2 , ...,x n )) € S' denotes 
the complex number transmitted by R and H x , H 2 ,..., H n 
respectively are the fading coefficients corresponding to the 
links R-Xi, R-X 2 , R-X n . The additive noises Zi, Z 2 „..., 
Z n are CAT (0, a 2 ). During the BC phase, R transmits a 
point from a signal set S given by a many to one map 
M Hl ' H2 '" ' Hn : S n — > S' chosen by R, depending on the 
values of H x , H 2 , H n . The cardinality of S' > 2^ n ~ 1 \ 
since X(n — 1) bits about the other (n — 1) users needs to be 
conveyed to each of X x , X 2 ,..., X n . 

A cluster is the set of elements in S n which are mapped 
to the same signal point in S by the map J\A Hl ' H2 ■■■■' Hn . 
LetC Hl H2 > - ' Hn = {Ci,C 2 , ...,Ci} denote the set of all such 
clusters. 

Definition 3: The cluster distance between a pair of clusters 
Li, Cj E c H i, H 2, --,H n ^ as gj ven m Q on me next page, is 

the minimum among all the distances calculated between the 



points (xi, x 2 , . 



€ C t and (x\,x 2 , 



t ) € Cj in the 



effective constellation seen at the relay node R. The minimum 
among all the cluster distances among all pairs of clusters of 
a clustering C Hl ' H2 ' - ' Hn is the minimum cluster distance of 
the clustering, as given in <|5j on the next page. 

During the MA phase, the performance depends on the 
minimum cluster distance, while during the BC phase, the 
performance is dependent on the minimum distance of the 
signal set S . Distance shortening, a phenomenon given in 



1 13 1, is described as the significant reduction in the value of 



d min (Hi, H 2 ,...,H n ) for values of (Hi, H 2 , .... H n ) in the 
neighborhood of the singular fade states. If the clustering used 
at the relay node R in the BC phase is chosen such that 



(C Hl 



M-2 



is non zero, then the effect of distance 
shortening can be avoided. 

A clustering C Hl < H2 '- - Hn is said to remove a singular fade 
state (H U H 2 , ...,H n ) e H, if d mm (('" " " ) > 0, i.e., 
any two message sequences (xi,x 2 , x n ) € S n that coincide 
in the effective constellation received at the relay during the 
MA phase is in the same cluster of C Hl ' H2 ' "' Hn . So, removing 
singular fade states for a n-way relay channel can alternatively 
be defined as: 



d m i n (Bi,H2,—,Hn) = min ^ ^ [Hi fan - + H2 (x2 

(xi ,X2 ,...,5Bn) iC^l 1^2 '■■■'^n)^ 4 ^' 1 

1 ) IC2 1 ■ • " J 71 ) ^ ('''I !^2' ,-,, '^ti^ 


- 1E2) + ••■ + -Hn (zn - a4j | 


(3) 


<filln 3 i H l> H a> #») = min 1^1 (^1 " ="0 + H 2 (*2 ~ 4 


) + ... + H n (x n -x' n )\ 


(4) 


dmin (C^^- -^") = min - + 

(x^ ^ ;l^2}■■'}3'7~t)}(3'•^_^■^•2'''" , "^' ^ ^ *^ } 
M H 1> H 2> ■••< H n ( Xl ,x 2 ,...,x n )^M H: L' B2 H " (k'j.X^,...,^) 


H 2 [x 2 ~ x' 2 ) + ... + H„ (x n - x' n ) | 


(5) 


dmin (Ci( firi ' Ha '-' H »5>,Hi ) H 2 ,...,fl'„) = min \H X (xi - si) 


+ H 2 (x 2 - x' 2 ) + ... + H„ (x n - x' n ) 


(6) 


M H l' H 2'"-' H ™(ai > X2,...,;z n ) ? 4A1 H l'- £f 2 


■■><) 




^ 1 ^ 2 ,...,^,.., JJ „ (xiiX2j ... i;Efci ... iXn) ^^ 1 ^,...^,...^ 




(V) 


where, x k = x' k , (x%, x 2 , ifc_i,a;fc + i, i„) ^ (4,a;2, ...,x' h _i,x' k 


x^J , Vii, 4 1 x 2i x' 2 , ...,x n ,x' n & S for k = 1, 2, .. 


., n. 



Definition 4: A clustering C Hl ' H2 '"' H " is said to re- 
move f/?e singular fade state (Hi, H 2 , H n ) € "H, if 
any two possibilities of the messages sent by the users 
(xi,x 2 , ■ £ n ) , (x\, x' 2 , zJJ € 5" that satisfy 



-ffi^i + H 2 x 2 



H n x n = H\x\ + H 2 x' 2 



H n x n 



■Hn 



, H, 



IS 



are placed together in the same cluster by the clustering. 

We denote the clustering which removes the singular 
fade state (Hi, H 2 , H n ) by C^ Hl ' H2 -- H ^ (selecting 
one randomly if there are multiple clusterings which re- 
move the same singular fade state (Hi, H 2 , ■ H n )). Let the 
set of all such clusterings be denoted by C-h, i.e., C-h = 
{Ci(^,-^)} : (Hi,H 2 ,...,H n )eH}. 

Definition 5: The minimum cluster distance of the clus- 
tering c{(-ffi>ff2,..,-ff„)} f or (Hi,H 2 ,...,H n ) e H, when the 
fade state (Hi, H 2 , H n ) occurs in the MA phase, denoted 
by d min (C {{Hl ' H ^-' Hn ^ ,H X ,H 2 , ...,H n ), is the minimum 
among all its cluster distances. 

If (Hi,H 2 ,...,H n ) £ H, the clustering C Hl ^ 2 
chosen to be C^ Hl ' Bl '—' Hn ^ G C u , that satisfies, 
d mm (C« H ^ H2 -> H «»,Hi,H 2 , H n ) > 

d min (ci( H ^' <)},H X ,H 2 ,. 
V(Hi,H 2 ,...,H n ) ^ (H[,H! 2 ,...,H' n ) e U. The clustering 
used by the relay is indicated to X\, X 2 ,...,X n using overhead 
bits. In order to ensure that Xk] k = 1, 2, .., n is able to decode 
the message sent by X\\ I = 1,2, ,,,n\ I ^ k, the clustering 
C should satisfy the exclusive law, as given in Q. We explain 
Exclusive Law in more detail in the next section. 
The contributions of this paper are as follows: 

• We propose a scheme that enables the exchange of 
information in the wireless n-way relaying scenario when 
A/-PSK is used at the n user nodes with totally two 
channel uses while attempting to remove the harmful 
effects of fading, extending the schemes given in (2|, J3], 
fTT) for n = 2, 3, 4 respectively. 

• For this scheme, the singular fade spaces are identified, 
enumerated and classified based on their removability in 
the given scenario. 

The remaining content is organized as follows: Section III 
demonstrates how a n-fold Latin Hyper-Cube of side M can 
be utilized to represent the network code that satisfies the 
exclusive law for n-way relaying when M-PSK is used at 



the end nodes. In Section IV we describe and enumerate the 
singular fade subspaces for the given scenario and in Section 
V, focus in on the removal of such singular fade subspaces 
using n-fold Latin Hyper-Cube of side M, Section VI provides 
some insights using simulations and Section VII concludes the 
paper. 

III. The Exclusive Law and Latin Hyper-Cubes 

The clustering C that represents the map used at the relay 
should satisfy the exclusive law [2| in order to ensure that 
Xk; k — 1,2, ...,n is able to decode the message sent by 
Xi; I = 1,2, ...,n; I 7^ k, where we assume that the nodes 
Xi, X 2 , X n transmit symbols from the A/-PSK constella- 
tion. Consider a M x M x ... x M (n times) array, containing 
M n entries indexed by (x\,x 2 , ...,x n ), i.e., the n symbols 
sent by Xi, X 2 , X n in the MA phase. For k = 1,2, n, 
fixing the k th dimension of this M x M x ... x M (n times) 
array, the M (n — 1) dimensional arrays obtained, denoted by 
say C l k , I = 1,2, M, are indexed by the M values taken 
by Xk- For fixed values of k and /, the repetition of a symbol 
in C k results in the failure of the k th exclusive law given by 
|7}. Thus, for the exclusive law to be satisfied, the cells of this 
array should be filled such that the M x M x ... x M (n times) 
array so obtained, is a n-fold Latin Hyper-Cube of side M, 
with entries from Z t = {0, 1, t - 1} for t > M^ n ~^ 
(Definition 1). The symbol Ci of a particular clustering 
{£i, ■■■,C t ] denotes the cluster obtained by putting together 
all the tuples (x\, x 2 , x n ) € S n such that the entry in 
the (xi, x 2 , x n )-fh slot is the same entry i from Z f . The 
adjoining figures Fig. [2] and Fig. [3] show the exclusive law 
condition for the three-way and four-way relaying scenario 
when 4-PSK is used at end nodes. 
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Xa=0 Xa=1 Xa=2 Xa=3 

Fig. 2. A 4-fold Latin Hyper-Cube of side 3 represents the exclusive law 
constraint for the relay map when 4-PSK is used at end nodes 

IV. Singular Fade Subspaces 

Definition 6: A set {(xi,x 2 , ...,x n )} € S n consisting of 
all the possibilities of (xi, x 2 , x n ) that must be placed 




Fig. 3. A 4-fold Latin Hyper-Cube of side 4 represents the exclusive law 
constraint for the relay map when 4-PSK is used at end nodes 



in the same cluster of the clustering used at relay node R 
in the BC phase in order to remove the singular fade state 
(Hi, H2, H n ), is referred to as a Singularity Removal 
Constraint for the singular fade state (Hi, H2, H n ) for n- 
way relaying scenario. 

At the end of the MA phase, the relay node receives a 
complex number, given by ([TJ. Using the ML estimate of 
this received complex number, R transmits a point from the 
constellation S' with cardinality at most M n , Instead of R 
transmitting a point from the M n point constellation resulting 
from all the possibilities of (xi, X2, x n ), depending on the 
fade states, the relay R can choose to group these possibilities 
into clusters represented by a smaller constellation, so that 
the minimum cluster distance is non-zero, as well as all the 
users receive the messages from the other (n — 1) users, i.e., 
the clustering satisfies the exclusive law. We provide one such 
clustering for the case of n-way relaying in the following. 

Suppose the fade coefficient in the MA phase, 
(Hi, H2, H n ), is a singular fade state, and T 
is a singularity removal constraint corresponding 
to the singular fade state (Hi,H 2 , ...,H n ). Then 
there exist (x\, x 2 , —, x n ), (x[, x' 2 , x' n ) € I\ 
(xi,x 2 , —,x n ) 7^ (x[,x 2 , x' n ) such that, 



Hixi + H 2 x 2 + ••■ + H n x n = Hix[ + H 2 x' 2 + ... + H n x' n 
>H 1 (x 1 - x[) + H 2 (x 2 - x' 2 ) + ... + H n (x n -<) = 



>(H u H 2 ,...,H n ) G 



x 2 — x' 



0) 



where for a n x 1 non-zero vector v over C, 

(v) 1 - = {w = (w 1 ,w 2 , w n ) e C" I toi^i + w 2 v 2 + ... + m„t)„ = 0} . 

(10) 

Note that wiUi + W2V2 + ■■■ + w n v n = J2i w i v i i s tne dot 
product over C (and not an inner product over C). 

Lemma 1: For a n x 1 non-zero vector v over C 
a (n — 1) -dimensional vector subspace of C n over 

Proof: Let w = (wi, W2, w n ) G (v) where v = 
(vi, i>2, v n ). Then, by definition, 



(v) is 



W1V1 + w 2 v 2 + ... + w n v n = 
=>wivi = -(w 2 v 2 + ... + w n v n ) 
^Wi = -Vi 1 (w 2 V2 + ... + w n v n ) 

=>(wi,W 2 , -,W n ) = (-Vi 1 (w 2 V 2 + ... + W n V n ),W2, ...,w n ) 
=>(W1,W2, ..;W n ) = w 2 (~v^ 1 v 2 , 1,0, ...,0) + 

^(-^r 1 ^, 0, 1, 0) + ... + w n (-Vi l v n , 0, 0, 1) 
(v) 1 - = span {(-Vi 1 v 2 , 1,0, ...,0), (-uf 1 u 3) 0, 1, ...,0), ... 

■ (-v^ 1 v n ,0,0, ...,1)} over C. 



Proving that the following subset of C" 

{ {~v- 1 v 2 , 1, 0, 0), (-v^v^, 0, 1, 0), (-v- 1 v n , 0, 0, 1)} 
is a linearly independent set over C of cardinality (n— 1) would 
be sufficient to prove that (v) ± is a (n— 1) -dimensional vector 
subspace of C". Let cti, a 2 , a n -i € C such that 

a 1 (-v~ 1 v 2 ,l,0, ...,0) + a 2 (-v- 1 v 3 ,0, 1 0) + ... 

... + a n _ 1 (-v- 1 v n ,0,0, ...,1) = (0, 0,..., 0) 
=>(— aiv~ 1 v 2 - a 2 v~ 1 vs — ... — a n v^ 1 t; n , a\, a 2 ,..., 

= (0, 0, 0). 



Comparing the 2 nd , 3 rd , n th components of the LHS and 
RHS, we get, ai = 0, a 2 = 0, a„_i = 0. 

■ 

Since xi, x 2 , ... , x n , x'i, x' 2 , x' n € S, where S is finite, 
there are only finitely many possibilities for the right-hand 
side of {9). Thus the uncountably infinite singular fade states 
(Hi, H2, H n ), are points in a finite number of (n — 1)- 
dimensional vector subspaces of C n over C. We shall refer to 
these finite number of vector subspaces as the Singular Fade 
Subspaces (3). 

We now give a detailed description of all the possibilities of 
singular fade subspaces for n-way relaying scenario when M- 
PSK is used at the end nodes. For the proof of the following 
Theorem, we extensively use the following Lemmas given with 
proofs in p4) . 

Lemma 2: Let AS denote the difference constellation of 
the signal set S, i.e., AS — {sj — s'Js^s^ 6 S}. Then, for 
any Af-PSK signal set, AS is of the form, 



A5 = {()) i.J j2sin(^)^| I is odd 



(11) 



U 



„ . , 7rZ 

2 sin — )e 



i( 5 



is even 



where 1 < I < M/2 and < k < M - 1. 

As a result of the above Lemma, the non-zero points in AS 
lie on M/2 circles of radius 2 sm(wl/M), 1 < I < M/2 with 
each circle containing M points. The phase angles of the M 
points on each circle is 2kn/M, if / is odd and 2kir/M+ir/M 
if I is even, where < k < M — 1. 



2 The proof is straightforward, yet given here for the sake of completeness. 



Lemma 3: [3] Let i\,i2,...,ih be the ordered indices 
corresponding to the non-zero components in Ax and Ax' (the 



location of non-zero components is the same in the vectors Ax 
and Ax'). For M-PSK signal set, \Axi\ = c|Ax-|,Vl < i < n, 
for some c e C, if and only if the magnitudes of the non-zero 
components in Ax are equal and the magnitudes of the non- 
zero components in Ax' are equal, i.e., lAxjJ = |Axj 2 | = 
... = |A^| and |Aa;JJ = |A<,| = ... = 

From Q and Lemma 2, the singular fade subspaces are 
given by, 

TTin 



XI 


- x > } - 

x 2 




X2 


H 


_ %n 


- x' 

- 





J»»2 




2sin ( 2 M ) e 



2 flinty" 

sin(^i)eJ m i 
sin(^)eJ m 2 



sin(^f)e^ 



where mi = 2fcj7r/M if is odd and 2kiii/M + ir/M if is 
even, where < hi < M — 1 for i = 1, 2, n. 

Theorem 1: There are £fe=i(5?) [(f) fc - (4f) + l] M fe_1 
Singular Fade Subspaces for n-way relaying when M-PSK 
constellation is used at the end nodes. 

Proof: The Singular Fade Subspaces are of the form 
([Axi, Ax2, •••) Ax n ]) where Ax/c <E AS, fc = 1,2, 
Let be the number of non-zero x^s. We fix the relative 
phase vector of the vector [Axi, Ax2, Ax n ] = w (say). 
The points in AS lie on M/2 circles. So there are (M/2) fc 
possibilities for absolute values of the non-zero components of 
w. There are M/2 possibilities for the case that the absolute 
values of all the components of w that are non-zero, are equal. 
From Lemma 3, the Singular Fade Subspaces resulting from 
all of these M/2 cases are the same, and hence account for 1 
out of the (M/2) fc cases. So for a fixed relative phase vector, 
there are [(M/2) k - M/2 + l] Singular Fade Subspaces. 
From Lemma 3, fixing the absolute values of the non-zero 
components of w, each distinct relative phase vector corre- 
sponds to a distinct Singular Fade Subspace. There are M fe_1 
distinct possibilities for the relative phase vector. So, there are 
[(M/2) fe - M/2 + l] M fe_1 Singular Fade Subspaces when 
w has k non-zero components. Here, k can take values from 
1 to n. Summing over all possible values of k, we have 
Efe=i(fe) [if) k ~ (f) + !] Mn ~ l Singular Fade Subspaces 
for n-way relaying when M-PSK constellation is used at the 
end nodes. ■ 

The above theorem coincides with the results given for 
n = 2 in fl4) , and the results obtained using explicit 
enumeration for n = 3 in |3| and for n = 4 in (TTJ. 
For illustration, we discuss the case when n = 5. In five- 
way relaying, user nodes (say) A, B, C, D and E trans- 
mit xa, xb, Xc, xd and xe € S respectively in the 
first channel use. Suppose the fade coefficient in the MA 
phase, (Ha, Hb, Hq, Hd, He), is a singular fade state. Then 
there exist (xa,xb,xc,x d ,x e ), (x' a , x' b , x' c , x' d , x' e ) G S 5 , 
(xa,xb,xc,xd,xe) ^ (x' A ,x' B ,x' c ,x' D ,x' E ) such that, 



H A x A + H B x B + H c x c + H D x D + H E x E 

= Hax' a + H B x' B + Hqx'q + H D x' D + H E x' E 
>H a {xa - x'a) + H B (x B - x' B ) + H c (xc - x' c ) 

+ H D (x D - x' D ) + H E (x E -x' E )=Q 
i 

xa — x , 



>(H A ,H B ,H C ,H D ,H E ) e 



xc ■ 
xd 
x E ■ 



F 



(12) 



The adaptive network coding for five-way relaying attempts 
at removing the singular fade subspaces for the case given by, 

xa ■ 



x B 
xc 
xd 
x E 



A 

$ 

1 D 



In the second channel use, relay node R 



transmits x r using a network coding map that depends on the 
values of x"a, x"b, xc, x"d and x"e- As explained in Section 
III, using a network coding map represented by a 5 -fold hyper 
latin-cube of side M (when M-PSK is used at the end nodes 
A, B, C, D and E) ensures that exclusive law is satisfied. It 
can be shown using explicit enumeration, that when the end 
nodes use 4-PSK, there are 13981 singular fade subspaces 
for five-way relaying, which coincides with Theorem 1, for 
n = 5, M = 4. 

V. Removing singular fade subspaces 
We cluster the possibilities of (x\, x%, x n ) into a cluster- 
ing that can be represented by an n-fold Latin Hyper-Cubes 
of side M, to obtain a clustering that removes the singular 
fade subspaces, and also attempts to minimize the size of 
the constellation used by R. This clustering is represented 
by a constellation given by S', which is utilized by the 
relay node R in the BC phase. This is done by first con- 
straining the M n possibilities of (xi, X2, x n ) transmitted 
at the MA phase, to remove the singular fade subspaces, and 
then using these constraints, filling the entries of an empty 
M x M x ... x M (n times) array representing the map to be 
used at the relay. This partially filled array is completed so as 
to form a n-fold Latin hyper-cube of side M, The mapping to 
be used at R can be obtained from the complete Latin hyper- 
cube keeping in mind the equivalence between the relay map 
that satisfies the exclusive law with the n-fold Latin Hyper- 
Cube of side M as shown in Section III. 

During the MA phase for the n-way relaying scenario, 
nodes X\, X2,---, X n transmit to the relay R. Let the 
fade state (Hi, H 2 , H n ) denote a point in one of the 
Efc=i(fe) [(fO " (t) + 1] M " _1 singular fade subspaces 
(Section IV). The constraints on the M n array repre- 
senting the map at the relay node R during BC phase 
for a singular fade state, can be obtained using the vec- 
tors of differences, viz., [xi — x\, x% — x' 2 ,---, x n — x' n ] 
contributing to this particular singular fade state. So, 



if (H u H 2 ,...,H n ) G ([X! 



Xi, X2 



u 2 , 



then, for (xi,X2, 
H 2 x 2 + ... + H n 3 



, x n ) , (x 1; x' 2 , x' n ) G S n , H\X\ + 



H lX [ + H 2 x' 2 



H n x' n . For a 



Algorithm 1: Obtaining the n-fold Latin Hyper-Cube of 
side M from the constrained M x M x ... x M (n — times) 
array 

Input: The constrained M X M X ... X M (n — times) 
array 

Output: A n-fold Latin Hyper-Cube of side M 

representing the clustering map at the relay 

1 Start with the constrained M x M x ... x M (n — times) 
array A" 

2 Initialize all empty cells of X to 

3 The *2j ---J *n)* ce H °f ^ i s me i*' 1 transmission of 
X\, the i| transmission of X 2 , the i„ transmission 
ofX„. 

for 1 < h < M do 
for 1 < i 2 < M do 



4 

5 

6 
7 
8 
9 

10 

11 

12 
13 
14 
15 
16 
17 

is end 



for 1 < i n < M do 

if cell (ii,«2; ■■■,««) of X is NULL then 
Initialize c=l 

if C c does not occur in the 

th 



(ii, i 2 , i n ) " cell of X then 
replace at cell (ii,i 2 , ■■■,i T . 
with C c ; 
else 
I c=c+l; 
end 



of X 



end 



end 



end 



clustering to remove the singular fade state (Hi,H 2 , ...,H n ), 
i.e., for the minimum distance of the clustering to be greater 
than (Section II), the pair (xi, £2, ■ ■■,x n ) , (x^, x' 2> x' n ) 
must be kept in the same cluster. Alternatively, we can say 
that the entry corresponding to (x\, X2, x n ) in the M n array 
must be the same as the entry corresponding to (x[, x 2 ,...,x' n ). 
Similarly, every other such pair in S n contributing to this 
same singular fade subspace must be kept in the same cluster. 
Apart from all such pairs in S n being kept in the same 
cluster of the clustering, in order to remove this particular 
fade state, there are no other constraints. Consider the ordering 
given as follows on the entries of the constrained M n array: 
(ci, c 2i on) < (c'j, c' 2 , c' n ) if Ci < where i is the first 
component among the n components, where c; ^ This 
constrained M n array can then be completed by simply filling 
the first empty cell in this order, with £j, i > 1 in the 
increasing order of i such that the completed array is an 71- 
fold Latin Hyper-Cube of side M (Algorithm 1). The above 
clustering scheme, however, cannot be utilized to remove all 
the singular fade subspaces, as shown in the following lemma. 

Lemma 4: The clustering map used at the relay node R 
cannot remove the singular fade spaces ([Axi, Ax 2 , Ax n ]) 
where at least one of Axk — for some k = 1,2, ...,n and 



simultaneously satisfy the mutually exclusive law. 

Proof: Let S — ([xi — x[, x 2 — x' 2 , x n — x' n }) ± be a 
singular fade state where for some 1 < k < n, Xk — x' k = 

0. Then, in order to remove S, (xx, x 2 , X)., x n ) and 
(x'n x' 2 , Xk, •••) x' n ) that must be kept in the same cluster. 
This would imply user not being able to distinguish 
between the messages xi and x\ for some 1 < I < n, I 7^ 
k, xi 7^ x\ sent by user Xi. This will clearly violate the 
mutually exclusive law, since in order to satisfy the mutually 
exclusive law, for the same value of Xk, all possible n-tuples 
of messages must be kept in different clusters. These two 
statements cannot be satisfied at the same time, hence such 
a singular fade subspace cannot be removed if the mutually 
exclusive law has to be satisfied by the relay map used in the 
BC phase. ■ 
We refer to the singular fade subspaces whose harmful effects 
cannot be removed by a proper choice of the clustering, as 
the non- removable singular fade subspaces also talked about 
in (14|. 

Corollary 1: There are [(f) n ~ (f) + l] M"" 1 
Removable and £fe=i(fc) [(f ) k - + l] M k ~ x Non- 
Removable Singular Fade Subspaces for ri-way relaying when 
A/-PSK constellation is used at the end nodes. 

Corollary 2: The number of non-removable singular sub- 
spaces is 0(Af I1_1 ) while the number of removable singular 
fade subspaces is 0(M n ). 

Thus, the number of non-removable singular fade subspaces 
is a small fraction of the total number of singular fade 
subspaces. For the five-way relaying scenario described in 
the previous section, there are 13981 singular fade subspaces 
for five-way relaying, out of which the scheme given in this 
paper removes 7936 singular fade subspaces using 5-fold Latin 
Hyper-Cubes of side 4. This can be done, as described above, 
by first marking the singularity removal constraints in the 
empty 4x4x4x4x4 array and then completing the array 
to form a 5-fold Latin Hyper-Cubes of side 4 using Algorithm 

1. We now illustrate the removal of a singular fade state for 
five-way relaying with the help of the following example. 

Example 1: Consider the case for five-way relaying sce- 
nario where 4-PSK is used at end nodes A, B, C, D and E, 
and one of the singular fade subspace to be removed is given 

by, 



s" = 



-1-3 




1-3 


-2j 




2 


-2j 


)■( 


2 


1 - j 




1 + 3 


1 + 3 




. - 1 + 3 . 





1+3 




-1 + 3 




2j 




-2 




2j 


H 


-2 




-1 + 3 




-1 - 3 




-1 - 3 




1-3 



The first vector is [-1-J, -2j, ~2j, l-j, 1 + j]. 
Now, can be obtained either as a difference of xa = — 1 

and x' A = j or as a difference of xa = — j and x' A = 1; — 2j 
can be obtained only as a difference of — j and j; 1—j can 
be obtained as a difference of xp — 1 and x' D = j or as a 
difference of xp = — j and x' D = 1; 1 + j can be obtained 
as a difference of xe = 1 and xe — —j or as a difference 
of xe = j and xe = —1- Thus, the entries corresponding to 



x c = 
x D = 










x c = 
= 1 










x c = 
□=0 = 2 










*o = 
^r> = 3 













r 4 


£l9 


^23 


£ 2 7 





£75 


£ 79 


c 83 


£ 8 7 





C-2 


£l39 


£l43 


£l47 





£l97 


^201 


^205 


£ 209 


1 


£s 


^20 


£ 24 


^28 


1 


£76 


r 80 


£ 84 


£88 


1 


£q 


£l40 


£ 144 


£l48 


1 


£l98 


C 202 


£ 206 


£210 


2 


£l7 


£21 


£ 26 


^29 


2 


£77 


£ 8 1 


£ 8 5 


£89 


2 


■£l37 


£141 


£■145 


^149 


2 


£l99 


C 203 


^207 


£211 


3 


£l8 


C 22 


£ 26 


£ 30 


3 


£78 


£ 8 2 


£86 


£90 


3 


r 138 


£l42 


£l46 


£l50 


3 


£ 200 


C 204 


£ 208 


£212 


x c = 1 
x D = 










x c = 1 
*D = 1 










33 = 1 
x D =2 










*o = 1 

x D = 3 













£ 3 1 


Cl8 


£37 


^41 





c 91 


£ 14 


C 97 


£l01 





£l51 


£ 155 


£ 159 


£l63 





■^213 


^217 


£■221 


£225 


1 


r 32 


£35 


£38 


£42 


1 


£■92 


£95 


c 98 


£•10.2 


1 


r 152 


r 156 


£l60 


£l64 


1 


^214 


£ 218 


£ 222 


£226 


2 


£ 3 3 


£l2 


£ 39 


£43 


2 


C 93 


£ 10 


£ 99 


r 103 
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£l53 


£l57 


£l61 
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C 219 
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x c = 2 
x D = 










x c =2 
X D = 1 










x c =2 
x D =2 










*o = 2 













^45 


£49 


£ 53 


%7 





C 1Q5 


£l09 


£ 113 


£ 117 





£ 167 


C 171 


^175 


^177 





■^229 


C 233 


r 237 


£239 


1 


£46 


£ 50 


£54 


£ 58 


1 




£ll0 


£ll4 


£ll8 


1 


£l68 


C 172 


^11 


£l78 


1 


£■2 30 


c 234 


£9 


£240 


2 


£47 


£51 


£55 


£ 59 


2 


£l07 


£lll 


£ll5 


r 119 


2 


£l69 


C 172 


^15 


^179 


2 


^231 


C 235 


£13 


£241 


3 


r 4S 
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£ 56 
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£ 112 


£ll6 


r 120 
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£l70 


£■17 A 


£l76 


£l80 


3 


£ 232 


c 236 


C 238 


£242 


a; c = 3 

x D = o 





I 
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3 


a; c = 3 
x D = 1 
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2 


3 







1 
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3 


*C = 3 
x D = 3 





1 


2 


3 





£&1 


£&5 


^69 


£73 





£12] 
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£ 129 


£l33 
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£l89 
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£ 243 


C 247 


C 251 


£255 


i 


£ 62 


£&6 


£70 


£74 


1 


£■122 


£l26 


£l30 


£l34 


1 


£l82 


£l86 


£l90 


£l94 


1 


■^244 


^248 


£252 
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£63 


£&7 


£71 


£l 
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£-131 


£l35 




£l83 
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£l91 


^195 
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£253 


£3 


3 
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£ 68 
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£ 5 


3 


£l24 


£ 128 
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£ 136 


3 


r l84 


£ 188 


£l92 


£ 196 


3 


^246 


£ 250 


C 254 


£7 



TABLE I 

Example 1: Latin-Hyper Cube representing the Relay Map where x e = 0, entries xc's and x d 's are as mentioned, xa's entries are 
along the rows and xfl's entries are along the columns of each 4x4 matrix for fixed values of x c and x d . 
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2 
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x D =2 





1 
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£104 
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£65 
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£224 
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£246 


£250 


1 


£lS5 


£l81 
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3 


£51 
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£59 


£l40 
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£235 
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£241 


£l0 


3 


£l77 


£257 


£ 260 
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TABLE II 

Example 1: Latin-Hyper Cube representing the Relay Map where xe = 1. entries xa's and x b 's are as mentioned, xc's entries are 
along the rows and x d 's entries are along the columns of each 4x4 matrix for fixed values of x a and x b . 



{(-1, -3, -j, 1, 1), U,3,3,3, -i)}. {(-1. -3, -3, -3, 1), 0,3,3, 1, -j)}> 
{(-1) ~3, -J, 1) j), (3,3,3,3, -1)}. {(-1, ~3, ~3, (j,3,3, 1. -1)}. 

{(-ji -J) -3, 1, 1), (1,3,3,3, -j)}. {(~3, ~j, ~3, ~3, 1), (1,3,3, 1, 
{(-3, -j, -j, 1J), (1,3,3,3, -1)}. {(-J, -j, ~j, -3,3), (1,3,3, 1, 
must lie in the same clustering representing the network coding 
map used at the relay node in the BC phase. Similarly the 
constraints resulting from the other three vectors above can 
be obtained. Replacing with 0,1,2,3, we get 

constraints in the form of entries that must be kept the same 
in the 4x4x4x4x4 array representing the clustering. 
For instance, corresponding to {(— 1,— j, — j, 1,1), (3,3,3,3, — j)}, 
the entries {(2, 3, 3, 0, 0), (1, 1, 1, 1, 3)} are kept in the same 
cluster £1 as shown in bold in Tables I, II, III, IV. The con- 
strained Hyper-Cube of side 4 is completed using Algorithm 
1 to form a 4-fold Latin Hyper-Cube of side 4. The completed 
Latin Hyper Cube is as shown in Tables I, II, III and IV, where 
the constraints are marked in bold. 

Similarly, a (removable) singular fade subspace can be 



removed by first constraining the array representing the relay 
map, and then completing the constrained array using the 
provided algorithm. 

VI. Simulation Results 

Simulation results presented in this section identify the cases 
where the proposed scheme outperforms the naive approach 
that uses the same map for all fade states and vice verse. 
Here, the channel states are distributed according to Rician 
distribution and channel variances equal to dB and the frame 
length is 256 bits. Fig. [4] compares the SNR vs bit-error-rate 
curves for three-way, four-way and five-way relaying scenario 
for (a) the adaptive network coding scheme presented in this 
paper with (b) the non-adaptive network coding using two 
channel uses, in which the same array is used by the relay as 
an encoder for all channel conditions. The details of three-way 
and four- way relaying can be found in 1 3 1 and JTTJ. The non- 
adaptive network coding for three-way, four-way and five-way 
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£43 


£33 


£77 




^96 


£l00 


£95 


£275 


1 


£l62 


£l65 


£l58 


£200 


1 


£219 


£262 


£227 


^266 


1 


£40 


£44 


£34 


£78 


1 


£99 


£101 


£274 


£276 


2 


£ 159 


£l63 


£ 197 


r 201 




£210 


£263 


£228 


c 267 


2 


£37 


£41 


£36 


£75 




£97 


£102 


£91 


£277 


3 


C 160 


£l64 


r 198 


C 202 


3 


£ 221 


C 264 


C 265 


C 268 


3 


£38 


£42 


£35 


£76 


3 


L 98 


£273 


£92 


C 278 



TABLE III 

Example 1: Latin-Hyper Cube representing the Relay Map where x e = 2, entries xa's and x b 's are as mentioned, xc's entries are 

ALONG THE ROWS AND X D 'S ENTRIES ARE ALONG THE COLUMNS OF EACH 4x4 MATRIX FOR FIXED VALUES OF Xa AND X B . 



x c = 
x D = 





1 


2 


3 


x c =0 
x D = 1 





1 


2 


3 


x c =0 
x D =2 





1 


2 


3 


x c =0 
x D = 3 





1 


2 


3 





£228 


£l96 


£235 


£232 





^252 


C 276 


^281 


£l84 





£l32 


£120 


£266 


£124 





£312 


£eo 


£&8 


£64 


1 


£241 


^195 


£236 


£ 231 


1 


^272 


C 275 


^282 


£283 


1 


£l31 


£119 


£ 267 


£123 


1 


£313 


£59 


£&7 


£63 


2 


^239 


^240 


C 250 


£24t 


2 


^189 


C 278 


£229 


£lS5 




£130 


£134 


£268 


£301 


2 


£72 


£70 


£l06 


£66 


3 


£227 


^251 


£249 


£234 


3 


£l79 


£277 


£l86 


£182 




£129 


£133 


£300 


£302 


3 


£73 


£69 


£105 


£65 


x c = 1 
x D = 




1 






X D = 1 




1 






x c = 1 

x D = 2 




1 






X C = 1 
x D = 3 




1 









£l92 


^242 


£256 


£l8S 





£l78 


£5 


£287 


£289 





£303 


£7 


£ll2 


£l08 





£314 


£56 


£52 


£318 


1 


£l91 


£253 


£255 


£l87 


1 


£l76 


£\ 


£l83 


£290 


1 


£304 


£3 


£lll 


£107 


1 


£315 


£55 


£51 


£l38 


2 


£l94 


£l80 


£233 


£24i 




£284 


£286 


•^181 


£l74 


2 


£ll8 


£114 


£l22 


£110 




£58 


£316 


£62 


£48 


3 


£l93 


£l90 


c 230 


£247 


3 


^285 


£175 


■^288 


£l73 


3 


£117 


£113 


£l21 


£109 


3 


£57 


£317 


C 61 


£47 


x c = 2 
x D = 





1 


2 


3 


x c =2 
x D = 1 





1 




3 


x c = 2 
x D =2 





1 


2 


3 


x c =2 
x D = 3 





1 


2 


3 





^254 


C 259 


C 262 


C 263\ 





£ 165 


£l62 


£■200 


£l58 





£102 


£99 


£307 


£309 





£44 


£40 


C 78 


£34 


1 


£257 


£212 


£261 


£ 265 


1 


^291 


£l61 


£l99 


£l57 


1 


£100 


£306 


£12 


£310 


1 


£43 


£39 


£l0 


£33 


2 


^218 


L 269 


C 264 


£ 273 




£l64 


£l60 


£s 


£l98 


2 


£305 


£98 


£l6 


£92 


2 


£42 


£38 


£l4 


£319 


3 


^258 


£219 


£219 


£274 


3 


£l63 


£159 


£4 


£l97 


3 


£101 


£97 


£308 


£92 


3 


£41 


£37 


£ 2 


£320 


x c = 3 
x D = 





1 


2 


3 


x c =3 
x D = 1 





1 




3 


x c =3 
x D =2 





1 


2 


3 


x c =3 
x D = 3 





1 


2 


3 





£211 


£207 


£279 


£ 204 





£ 150 


c 294 


£l48 


£298 





£90 


C 86 


£94 


£83 





£30 


£26 


£22 


£115 


1 


^208 


C 166 


C 269 


£ 203 


1 


£l49 


c 295 


£l47 


£299 


1 


£39 


£13 


£93 


£si 


1 


£29 


£25 


£21 


£l5 


2 


^209 


£224 


£214 


£ 279 


2 


£292 


£ 296 


£l56 


£l42 


2 


£ss 


£84 


£so 


£311 


2 


£28 


£321 


£32 


£18 


3 


C 296 


£295 


£271 


£ 280 


3 


£293 


£ 297 


£l55 


£l41 


3 


£37 


C 83 


£79 


£9 


3 


£27 


£322 


£31 


£ll 



TABLE IV 

Example 1: Latin-Hyper Cube representing the Relay Map where x e = 3, entries xa's and x b 's are as mentioned, xc's, entries are 
along the rows and x d 's entries are along the columns of each 4x4 matrix for fixed values of x a and x b . 

relaying utilizes the same 3-fold Latin Hyper-Cubes of side 
4-fold Latin Hyper-Cube of side 4 and 5-fold Latin Hyp 
Cube of side 4 (respectively) for all channel conditions. 

For any communication system, the effects of additive no 
are predominant at low SNR and the effects of noise due 
fading dominate at high SNR. Since adaptive network codi 
scheme attempts at reducing the effects of fading for the 
way relaying scenario by removing a part of singular ft 
states, (a) performs better than (b) at high SNR. This a 
implies that the performance of the scheme in the MA ph; 
(removal of singular fade spaces takes place in the MA pha 
is predominant at higher SNRs, and the performance of i 
scheme in the BC phase dominates at lower SNRs. Also, 
the number of user nodes n increases, the SNR at which i 
performance of adaptive network coding improves over t 
performance of non-adaptive network coding increases, as c 
be seen in the plot, since the size of the received constellati 
in the MA phase and hence the size of the constellation us 



:he 
:an 
on 
;ed 




m SiT-Ai^iu 'n-n I U;ir i'ic:r£U-ri£ r*" ..lie el 

fir tive-w:r, rd:miu ivhc I 1-PS\ i- n-.ee :<l ;ne I wits- 
# Adaptive PIiym.1 Layer \e v,vil (Vehie using too diannel 

lb: livE-wi\ lelayiriJ ivNen 4-l'S.< i- n-cc ill ;m: uiiles 
. Non-Adaptive Fhy> a\ l.a.wr Nclworl; (liuliri* using too channel i 
for fonr-way rdj\ ing 4-PSK is .isd al end nodes 
Adaptive Physical Layer V wml del ne nsicig two channel uses 
ib tcur-e;ay id i.., . I ■■ PSKi iscd al nid nodes. 
j^Noti- Adaptive Physical Layer Network Coding using too channel 

for three m\ 1 1 ised at endnodes 

_^_Adaptive Physical Layer Network Codingnsing two channel uses 
for three-way relaying when 4-PSK is used at end nodes 



Eb/No, dB 

Fig. 4. SNR vs BER curves for different schemes for multi-way relaying 
when the Rician Factors is 20 dB 



in the BC phase increases with increasing values of n. 



VII. Conclusion 

We consider the n-way wireless relaying scenario, where 
n nodes operate in half-duplex mode and transmit points 
from the same M-PSK constellation. Information exchange 
is made possible using just two channels uses, unlike the 
existing work done for the case, to the best of our knowledge. 
The Relay node clusters the M n possible transmitted tuples 
(x\, X2, x n ) into various clusters depending on the fade 
states such that the exclusive law is satisfied and some of 
the singular fade subspaces are removed. This necessary 
requirement of satisfying the exclusive law is shown to be 
the same as the clustering being represented by a n-fold Latin 
Hyper-Cube of side M. The size of the clustering utilizing 
modified clustering may not be the best that can be achieved, 
and it might be possible to fill the array with lesser symbols. 
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